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Abstract 

Four-photon scattering in nonlinear waveguides is an important physical process that allows 
photon-pair generation in well defined guided modes, with high rate and reasonably low noise. 
Most of the experiments to date used the scalar four-photon scattering process in which the pump 
photons and the scattered photons have the same polarization. In birefringent waveguides, vec- 
torial four-photon scattering are also allowed: these vectorial scattering processes involve photons 
with different polarizations. In this article, the theory of four-photon scattering in nonlinear, bire- 
fringent, and dispersive fibers is developed in the framework of the quantum theory of light. The 
work focusses on the spectral properties and quantum correlations (including entanglement) of 
photon-pairs generated in high-birefringence and low-birefringence fibers. 
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I. INTRODUCTION 



When two intense monochromatic pump beams (or quasi-monochromatic pulses) are 
launched together in an one-dimensional nonlinear medium, matter-light interaction 
results in the scattering of pump photons to other wavelengths. The main scattering process 
consists in the spontaneous conversion of two pump photons with angular frequencies uuoi 
and ujq2 into a red-shifted Stokes photon (signal) and a blue-shifted anti-Stokes photon (idler) 
with angular frequencies u s and u a , satisfying u i + ojq2 = ^>s + u a . This elastic process in 
known as a four-photon scattering (FPS). It conserves the number of particles as well as the 
total energy in the field. 

During the recent years, FPS has attracted much attention from experimentalists because 
it allows to generate correlated, and sometimes entangled, photon pairs in optical fibers 
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19 . l20l . l2l| . These fiber-optics 
photon-pair sources constitute a major progress in the development of quantum photonics 
because they offer the advantage of being compact guided-wave sources that can be easily 
connected (with only marginal losses) to standard transmission fibers. This constitutes 
an important advantage for applications in fiber quantum communication (quantum key 
distribution [22J, remote coin tossin 



FPS in silicon waveguides 
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23]) and fiber quantum computation [24J. Recently, 
271 ] has also been demonstrated and used to produce 
photon-pairs suitable for silicon-on-insulator quantum circuits. 

In contrast with the considerable amount of experiments, theoretical work on FPS has 



been rather limited 
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28l . l29l . |3(| . The early work of Wang et al. [28( gives a quantum 
description of FPS in the scalar case, where the pump photons and the photon pairs share the 
same polarization. That treatment is restricted to a single monochromatic pump (degenerate 
case, co>oi — ^02) an d was generalized later Uj to account for the spectral shape of short 



pump pulses. In addition, Ref. 
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introduced, for the first time, first-order perturbation 



theory for studying FPS. Perturbation theory turns out to be a very suitable approach and 



will be used extensivelly in this work. In 2006, Lin et al. |29l . 1301 ] examined the theory of 
Raman noise in fiber-optics photon-pair sources and observed that the impact of Raman 
noise can be reduced by producing photon pairs in birefringent fibers through a vectorial 
FPS which generates photons that are orthogonally polarized with respect to the pump. So 
far, vectorial FPS processes have not been investigated experimentally (except in [9| in a 
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reverse degenerate FPS configuration). 

In this paper, the theory of FPS in a nonlinear, birefringent, and dispersive fiber is 
developed in the framework of the quantum theory of light. The quantum state of light at the 
output of the fiber in computed using first-order perturbation theory. Only degenerate FPS 
from a single monochromatic pump (0% = cj 2 = ^0) is considered and Raman scattering is 
not included. The work focus on the spectral properties of the generated photon-pairs and 
the kind of correlations that can be obtained in vectorial FPS processes. 

If the optical power of the pump and the propagation length are large enough, the pho- 
tons that have been generated by FPS may be amplified by stimulated four-wave mixing 
(FWM). FPS and stimulated FWM are two different aspects of the same phenomenon. To 
highlight this point, the properties of vectorial FPS are derived from the quantum nonlinear 
Schrodinger equations (see Sec. lIIII) . in close analogy with the classical analysis of modulation 
instabilities (self phase-mached FWM processes in fibers). The cases of high-birefringence 
and low-birefringence fibers are treated separately (Sec. HV]and Sec. |V]), as in the classical 



theory 3jj. Although, FPS and classical FWM are strongly related, they are still physically 
distinct processes because, as will be shown in Sec. I VI I parametric gain is not required for 
FPS. Therefore, the spectrum of photons generated by FPS is usually wider than the range 
of frequencies experiencing parametric amplification due to classical FWM. 

II. QUANTUM THEORY OF NONLINEAR PROPAGATION 

Single mode optical fibers are one-dimensional propagation media that are weakly dis- 
persive, weakly nonlinear and possibly birefringent. Only linear birefringence is considered 
here. The propagating electric field operator (in the Heisenberg picture) can be written as 

E= E 3 (z,t) ej F(x,y), (1) 

where T{x,y) (ff \J r (x,y)\ 2 dxdy = 1) is a function of the transverse coordinates describing 
the profile of the guided mode (it is assumed to be the same for both polarizations and all 
the frequencies involved), e,,- (j e {x,y}) are unit vectors along the optical axes, and 
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Here kj(u) and nj(u) are the propagation constant and refraction index on the j-axis. 
The operators bj(z,u) are ordinary photon annihilation operators. In absence of any non- 
linearity, the operators bj(z,u) are ^-independent. However, the propagation of a strong 
monochromatic pump in the fiber (with angular frequency Uq) makes the different bj(z,u) 
^-dependant. For any u> ^ u) , the z-dependance can be split in two parts: (i) a FWM part 
that, in the quantum theory, also gives rise to the FPS phenomenon and (ii) a nonlinear 
phase shift due to the cross phase modulation (CPM) induced by the pump. To make the 
CPM effect explicit, we write bj(z,w) = b'Az,w) exp [i 4>cpm( z )], where 



<j )CPU (z) = 2 1 [P 0j + -P - j j. (3) 

The z-evolution ofb'j(z,uj) is only due to wave-mixing. The phase factor due to the CPM 
depends on the polarization of the pump. P j and P j are t ne pump powers propagating on 
the j-axis and the j'-axis, respectively (j = y if j = x and vice versa). The factor 7 is the 
non-linearity parameter of the fiber, defined in Sec. II III [Eq. ( fl8i) ]. 

Because the fields that will be considered have a small bandwidth (usually up to sev- 
eral tens of THz) around the central pump frequency uq, one can neglect the frequency 
dependence of the square root in the right-hand side of Eq. (T5]) and write 



Ej (z,t) = ,\—^- e-«** / aj (z,Q) = dn, (4) 

ze nj c J^ee ^271 



where n j0 = nj(u ) and Q = uj — u is the frequency detuning. We also set a,j(z,Q,) = 
b'Az,uj) and (3j{Q) = kj(u). The operators aj(z,Q) and aUz,Q) satisfy the usual bosonic 
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commutation relations: 



a j (z,n),a],(z,n')\ = 6 jjr 6(n-tf), (5a) 

a^(2,n), 0^(^,0')] = 0, (5b) 
aUz,n),aUz,n')] = 0. (5c) 



Note that the fiber is not infinite: it begins at z = and end at z = L. Since the linear 
dispersion relations are discontinuous at these two points, a mode-coupling between forward 
and backward propagating photons occurs (Fresnel reflection due to impedance mismatch). 
We will neglect these effects and consider that a,j(0,Q) = lim^-oo aj(z, Q) and aj(L,Q) = 
lim^oo dj(z, fl) (there is no non-linearity outside the fiber). 



In order to solve the nonlinear propagation problem, we need the unitary operator U(L, 0) 
that maps a,j(0,Q) on aj(L,Q) for any detuning Q: 

aj (L, Q) = U ] (L, 0) aj (0, Q)U(L, 0) . (6) 

Obviously, the evolution of aj from z = to z = L is a continuous process. Therefore, there 
exists an unitary evolution operator U(z 2 ,zi) such that 

a,j(z 2 ,Si) = U\z 2 , z 1 )a j (z 1 ,Q)U(z 2 , z x ) (7) 

for any couple of coordinates (zi,z 2 ) in [0, L]. The continuity in the variables z\ and z 2 

implies the existence of an Hermitian operator G(z) (the "infinitesimal generator") such 

that, for small 5z, U(z + 8z, z) = 1 + lG(z)Sz + o(8z). Using this last relation and Eq. (j7j), 

one finds that the annihilation operators dj(z, Q) satisfies the following Heisenberg equations: 

d 

ih— aj (z, Q) = [G(z),aj(z, Q)] . (8) 

As it will be shown later, Eqs. (jHJ) are just the usual coupled-mode equations of nonlinear 
optics. In principle, solving these equations yields the solutions a,j(L,Q). The operators 
aj(L, fl) and the initial quantum state of light provide a complete knowledge about the 
final state since the mean number of photons (il)\a,j(L,Q)aj(L,d)\ip) in each mode can be 
computed, as well as any photon correlation between modes. 

An alternative approach consists in working in the interaction picture instead of the usual 
Heisenberg picture. The passage from one picture to the other one is performed according 
to the transformation 

\^(z)) = U(z,0M, (9a) 
A\z) = U{z, 0)^)^ {z, 0) (9b) 

where is an arbitrary initial quantum state of light and A an arbitrary quantum operator 
in the Heisenberg picture. For instance, it follows from Eq. (|7j) that the interaction-picture 
annihilation operators are z-independant: 

afa, tt) = U(z, 0)aj{z, tt)U f (z, 0) = ^(0, Q) = ^(O). (10) 

As a general rule, the interaction-picture evolution of operators is due to dispersion, birefrin- 
gence (the linear properties of the fiber) and CPM (but not FWM), while the interaction- 
picture evolution of the quantum states is only due to the FWM: 

-ih^\tf{z)) = G?{z)\tf{z)), (11) 
5 



where 

G\z) = U{z, 0)G(z)rf(z, 0). (12) 

From (EH), we see that \ip) = |?/> 7 (z = 0)). Inserting (J9fei) into fTTTT) gives an equation for 
U(z,0): 

l-U{z^)= l -G I {z) U(z,0), 
which can also be written as an integral equation: 

U{L, 0) = 1 + l - [ L G\z') U{z', 0) dz'. (13) 
" Jo 

The main advantage of working in the interaction picture is that Eq. ( {TBI) can be solved 



iteratively using the standard Dyson's perturbation expansion 32| from quantum theory of 
scattering. For vanishing non-linearity, the second term of the right-hand side of (jTBl) can 
be neglected, so the zero-order approximation of U(L, 0) is just U^(L, 0) = 1. For weak 
non-linearity, a better approximation is obtained by replacing U(z',0) under the integral 
sign in Eq. ( jTBl) by the zero-order approximation U^°'(z',0) = 1. This yields the first-order 
approximation 

C/W(L,0) = 1 + ^ / G T (z r ) dz. (14) 
Jo 

Continuing the recurrence, better approximations can be found using: 

U {n \L, 0) = 1 + l - [ L G'iz') U^iz', 0) dz'. (15) 
" Jo 

Since a single coherent monochromatic pump field is launched in the optical fiber, the 
initial state \ip (z = 0)) is the vacuum state for all modes except the pump mode. In the 
article, we call this state |0). As we will show later, using Eq. (J9K) with the first-order 
approximation of the evolution operator U^(L, 0) yields the final state 

POO 

\tf(L)) = \0) + / U(L,n) |l£X n )dQ 
Jo 

POO 

+ / £ yy (L,n) iiJi.i^dn 

Jo 

poo 

Jo 

POO 

+ / ty.(L,n) |l»,l-n>di2, (16) 
Jo 

where |1^, l^ n ) = a](J7)at,(— fi)|0). Two photons are created as a consequence of the an- 
nihilation of two pump photons. This process is called a FPS. The two-photon amplitudes 
£jf(z, fl) determine the polarization and spectral properties of the created photon pair. 



Since a continuous-mode formalism is used, the states l^ n ) are not dimensionless 
(they have the dimension of time) and cannot be normalized in the usual sense. (Only 
linear combinations of a continuous set of these states are physical.) When the pump has a 
finite duration T (but is still long enough to allow the monochromatic-pump approximation), 
it is often easier to describe the scattering in terms of discrete modes, the spectral width of 
which is equal to Aa> = 27r/T. The transition from the continuous- mode to the discrete- mode 
description requires the substitutions 



dj(Q) VAlj — > aj(fl). 

In the discrete-mode description the two-photon states |1q, l^_ n ) are dimensionless and nor- 
malized. 

It is important to note that the whole perturbation scheme strongly depends on the way 
the operators aj(L, Q) are defined by Eq. (jlj). In particular, the fact that we factored out the 
CPM phase factor makes the perturbation scheme used in this work different from the one 
used in the previous perturbation-theory description of FPS Hj. The perturbation scheme 
used in this work gives more precise results than the approach in H[ at any perturbation 



order 



33]. 



Using the theory presented in the next sections, we will be able to compute the two-photon 
amplitudes £jji(z,Q) and show the connection between FPS and classical FWM. 



III. THE QUANTUM NONLINEAR SCHRODINGER EQUATIONS 

If the non-linearity is neglected, the fiber exhibits two linearly polarized optical modes 
with propagations constants (3 X and (3 y . Because of the dispersion, both propagation con- 
stants vary slowly with the angular frequency of light. In many cases, a second order Taylor 
approximation around the pump frequency uj is precise enough to capture the essence of 
the physics involved. We can write 

/3 J (fi)=^+A^ + ^ 2 , (17) 

with (3 n j = d n Pj/dQ n (Q = 0) and j G {x, y}. The two polarization modes have thus 
different phase and group velocities. The parameters A/3 = (3 Qx — j3 0y and Ap 1 = fi lx — j3i y 
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measure the phase and group velocity mismatch, respectively. We will assume that the 
group-velocity dispersion (GVD) parameter [p 2 j in Eq. (IT7|) ] is the same for both axes and 
write p 2 = 02x — @2y When 2tic/u>o is close to the so-called zero-dispersion wavelength, j3 2 
can be so small that higher-order terms in the Taylor expansion (|l7p must to be taken into 
account for consistency. Although this situation is encountered in many experiments that 
generate far detuned photon pairs using scalar FPS, we will not consider that case in this 
study of vector FPS. 

The intrinsic non-linearity of silica fibres is of third order and can be considered isotropic 



31] . As a consequence, the nonlinear susceptibility tensor has only one independent element 



Xxxxx- The non-linearity parameter of the fibre is defined by: 

where n and A e s are, respectively, the effective linear index and the effective area of the 
fibre at the frequency uq. The constants eo and c are the vacuum permittivity and velocity 
of light. 

Light propagation in a single-mode biref ring ent fiber is properly described by the a set 



of coupled nonlinear Schrodinger equations 31| for the slowly varying envelopes A x and A 



of the two polarization modes. These are connected to the electric field components by the 
relations 

Ej(z, t) = AT Aj(z, t) e^*-" *). (19) 

Fast time oscillation at the pump frequency ojq and fast space oscillations with the pump 
propagation constants (3q x and /3 0y have been factored out. In quantum theory, the enve- 
lope fields are operators. The dimensional constant Af = 1/ ^2e nj c is chosen so that 
(Aj(zo,t)AAzo,t)) is the mean optical power flowing through the plane z = z Q at time t 
with a polarization along the j-axis, j G {x, y}. If the approximations explained above hold, 
the coupled quantum nonlinear Schrodinger equations are 

+ A- d ~w + 4 1 = * 7 { A * A - + 1 A * A y) A ' + l l A * A < ^ 2iA *M 2to ) 

+ i3i » 8 ~w + ' § ^ = 1 7 ( A * A > + 1 • 4 '- 4 -) A « + ' I AlA > e+2 " aft,7 < 20b ) 

In the following, we want to describe how photons are scattered from an intense monochro- 
matic pump wave. The pump wave must be a known steady-state solution of Eq. fl20|) . In 



S 



principle, there is no need for the pump being polarized along an optical axis. Its polariza- 
tion state can rotate. The general expression of the steady-state solution of (1201) is known 



31] . However, it involves elliptic functions and is difficult to handel for the purpose of the 



present work. We will therefore restrict the analysis to particular cases of special interest. 
We will distinguish high-birefringence (HB) and low-birefringence (LB) fibers. In the first 
case, the pump can have any polarization state; in the second case, we will restrict the 
analysis to linearly polarized pump wave parallel to an optical axis to make the problem 
tractable. 



IV. HIGH-BIREFRINGENCE FIBERS 

In the HB case, the beat length 27r/A/3o is much smaller than any other relevant length 
scale: the phase factors in the last terms of the right-hand side of Eqs. ( 120]) oscillate so quickly 
that they simply average to zero. To a good approximation, Eqs. ( 1201) can be replaced by 

dA x 
dz 

dAy_ 

dz 

An initial classical monochromatic wave has a constant complex envelope 

A x (0,t) = y/Ptee* ', (22a) 
Ay(0,t) = v^e** (22b) 

Pq x and Po y represent the power of the optical fields polarized along the x- and y-axis, 
respectively. Eqs. ( 1211) admit a classical monochromatic solution that is compatible with the 
initial conditions (1221) : 

A x (z,t) = e ie °* e i7 ( Pte+ ^) 2 , (23a) 

A y (z,t) = y/P^ e i6 °y e^+i^K (23b) 

However, this solution is never stable: during the propagation, some pump photons are 
destroyed and new photons, at different frequencies, are created, making the field polychro- 

9 



-P: 



lx~ 



dA x & a 2 a, 

dt 1 2 dt 2 



+27 \AlA x + -A\A y 
d_A^_ (3 2 d 2 A y 

Ply ^ . 4 



dt 



2 dt 2 



r ( A\A y + -AlA 



A, 



Ay. 



(21a) 



(21b) 



matic. In order to demonstrate this point, let's introduce the ansatz 

A x (z,t) = (>//£ + «»(*,*)) e i9 °* J< p °* + I p °y> , 
A y {z,t) = (^/P^+u y {z,t)) e ie °y e K^<+§ p °*> 



(24a) 
(24b) 

in Eqs. (12~T|) . The fields Uj(z,t), j G {x, ?/}, can be seen as perturbations to the stationary 
solution f[2"31 . By injecting the ansatz (|24|) into Eqs. f[2~Tj) and retaining only the terms linear 
in Uj(z,t), we obtain the following equations: 

du x du x .(3 2 d 2 u x , n 



O ()l 2 +ilP0y(u y + ul) 



3 



dUy .^O^Uy 



Of 



— I- 



(25a) 



(25b) 



A. Coupled-mode equations 



We now make use of Eqs. (j3J), (J4|), (fT9|) . and (|24|) to write the mode expansion of the 
fields Uf 



Uj(z,t) 



~2T 



-i6 0j 



Inserting the expansion (1261) into Eqs. (1251) . we obtain the coupled-mode equations 



(26) 



d_ 

dz 



9^ 



3 



+ e' 



+ *g7V-P0X-P()l/ 



e i(e 0H -e .) e -i[-A/3in+ 7 (Poy-P 0a; )]z a ( Z; 0) 



+ giCeo^+fioj,) e -i[-A/3if7+/3 2 n 2 +7(Po I +Poj / )] 2 n t^ _( 



al(z,-Q) 



(27a) 



(27b) 



These equations show that the non-linearity couples the x- and y-polarized modes of fre- 
quency uj ± Q. Note that the evolution of a x (z, —fl) and a y (z, — fl) is obtained from fl2~TK) 
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and (127b ) by replacing Q by —Q. Comparing Eqs. (l2"Tj) with the general formula dHJ shows 
that in the HB limit, the infinitesimal generator is given by 



G hh (z) = ^ 7 P 0x ( e 2i0 °* J™ e- i ^ Q2+2 ^> at(z,Q)al(z,-Q) dQ + H.c.) 
+ ^P 0y (e 2ie °y J°° e -W2&+-yPo v )z al(z,Q)al(z,-Q) dQ + H.c.) 



+ fi-^ly/PoxPoy 



'OO 

e i(0a*-0oy) I e -i[A^n +1 (P 0x -P 0y )]z fl t ^ n) dfi + H.c. 



— oo 



2 / 

POO 

e i(0 Oy -e Ox ) / e -i[-A^n+^p 0y -p 0x )]z w ty ax ( Zj q) dfi + h.c. 



ft , 



'OO 

e *(^+M / e -i[^+p 2 ^{P 0x +p 0y )] z a f (Z) n)Q t _q) dfi + h.c 



3C 



^o,+M y°° e-il-AAn+ftn^+^+Po,)], a t(^) a t( Z; _q) dft + H.c.) (28) 



B. Four-photon scattering 



So far, we used the Heisenberg picture to derive the coupled-mode equations f[2T|) and the 
corresponding infinitesimal generator (T2"8"j) . However, as explained in Sec. [TTJ more physical 
insight is gained by working in the interaction picture, especially when Dyson's perturbation 
series technique is used to compute the quantum state of light at the output of the fiber. 
To apply this technique, we first need to find G^ b (z), the interaction-picture expression 
of the infinitesimal generator. Using Eqs. ( fT2l) and ([7j), one sees that G^ h (z) is obtained 
from Ghb(^) [Eq. 0281) ] by replacing all the annihilation operators aj(z, ±Q) by their values 
^•(0, ±Q) = cij(±n) at z = 0. Note that G[ h (z) ^ G hb (0). 

We now compute the quantum state of light at the output of the fiber using the first 
order perturbation theory: \ip(L)) = C/ (1) (L,0)|0), with U^(L, 0) given by Eq. flU). The 
result f|T6|) anticipated in Sec. [TT] is readily obtained, with 
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£,xx(L, fl) 
£yy(L, fl) 
£, xy (L, fl) 



i ( 7 P x^) e 2ie °> 



e -i(/3 2 n 2 +2 7 P 0l )f ginc 



i (jPoyL) e 



2i8 0y -i((3 2 n 2 +2'yPoy) 



2 sine 



(p 2 fi 2 + 2 7 p 0x .) 

{(3 2 fl 2 + 2 1 P 0y ) 



(29a) 
(29b) 



do x +eoy)„-i{Af3 1 n+f3 2 n 2 +'y(Pox+Poy)] 



vn-% 



x smc 



(Af3 1 fl + (3 2 fl 2 + 1 (P 0x + P 0y )) 



£ yx (L,fl) = i (jpy/PoxPoyLj e 



(29c) 



i(e Ox +0o y ) e i [A/3i Q-/3 2 Q 2 - 7 (P x +Po v )] f 



x smc 



(Aftfi - /3 2 fi 2 - 7 (P x + Pfly)) 2 



(29d) 



Eqs. f fT6|) and (129]) describe the FPS in high-birefringence optical fibers. Four different 
processes can be distinguished according to the polarization of the generated photons. They 
are summarized in Tab. [J together with the corresponding two-photon amplitudes £jj'(L, fl). 
Note that in the processes corresponding to the amplitudes £ XX (L, fl) and £ yy (L, fl) the pump 
photons are co-polarized and have the same polarization as the scattered photons (scalar 
scattering). In the processes corresponding to the amplitudes £ xy (L,fl) and £, yx (L,fl) the 
pump photons have orthogonal polarizations (vectorial scattering). Let's analyze the scalar 
scattering first. 





Stokes on x 


Stokes on y 


anti-Stokes on x 


fl) 


£ xy (L, fl) 


anti-Stokes on y 


£,yx(L, fl) 


£yy(L, fl) 



TABLE I: Possible scattering processes in high-birefringence fibers and the corresponding two- 
photon amplitudes 



C. Scalar scattering 

If the pump wave is linearly polarized along an optical axis of the fiber (the x-axis), 
only the scalar scattering corresponding to that axis takes place (£ yy (L,fl) = £ xy (L,fl) = 
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\tf(L)) = |o> + / n) i-n) da (30) 

This formula shows that the generated pairs are in an energy- entangled state: the energies 
of the signal and idler photons are correlated while, at the same time, each photon is in a 
coherent superposition of a continuum of possible energy eigenstates. 

Let's compute the average photon-flux spectral density f x (L,Q) at the output of the fiber 



for an arbitrary angular frequency ujq + Q. This is given by 34j 



f x (L, Q) — — lim - / / (V J (L)|4(fi 1 )a :r (fi 2 )|V / (L)}dn 1 dfi 2 . (31) 

27T £^0 £ _/ n _e 



Using Eqs. flUD and (TJSi), we find that 



(/3 2 ft 2 + 2 7 P 0;E ))^ 



(32) 



2tt 2tt 

The flux /^(-L, Q) is plotted in Fig. [1] as a function of the Q for a fixed pump power and 
different fiber lengths. Panel (a) and (b) are for normal and anomalous dispersion, respec- 
tively. Note that f x (L,Q) = f x (L, — O) because signal and idler photons are always created 
in pairs. The agreement between the first-order approximation (I32p [dashed black lines] and 
the exact solution [solid gray line, see Eq. (|57[) ] is excellent as long as the fiber length is 
much shorter than the non-linearity length: 

L « Z% = -i-, (33) 

i.e. f x (L, Q) -C 1 for any value of Q. In this limit, the sign of the dispersion influences only 
slightly the spectral shape of the created photons and the spectral width of the fluorescence 
spectrum (full first-zero width of the sinc-function) is Af2 sca i m 2a/27t/(|/5 2 |L). The fact 
that the sign of /? 2 has almost no influence on the FPS spectrum strongly contrasts with 
the conditions under which scalar modulation instability (the stimulated counterpart of 
scalar FPS) can be observed. The developpement of scalar modulation instability requires 
parametric gain, which only exists in the anomalous dispersion regime (if the dispersion 
relation is quadratic, as has been assumed through this work, see Eq. (TT71)). FPS and the 
stimulated FWM phenomena called modulation instabilities will be further compared in 
section PVTl 

If the pump wave has a finite duration T, the frequency space can be divide in modes 
of finite spectral width Au; = 2-k/T (see the end of Sec. [TTJ) . Formula fl32|) then shows that 
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FIG. 1: Photon-flux spectral density f x (L,Q) for the anomalous (a) and normal (b) dispersion 
regimes. The first-order approximation f x (L,Q) = \£ XX (L,Q)\ 2 / (2tt) (dashed black lines) is com- 
pared with the exact value (solid gray lines) obtained by solving Eqs. ([27|) in the Heisenberg picture. 
The parameters used in this figure are 7 = 3 /(W km), Pq x = 300 mW, $2 = ±20 ps 2 /km [— in 
panel (a) and + in panel (b)], and L takes successively the values 100 m, 200 m, and 300 m. 

\£ XX (L,Q)\ 2 is the mean number of photons n x (L,Q) in the mode ujq ± fl If (1331) holds, 
the mean number of photons in any given mode u + Q is much lower than one and is 
therefore numerically equal to the probability of finding a photon in that mode. The limit 
lim^oo n x (L, £l)/T = p xx (L, Q)dQ = \l; xx (L, Q) \ 2 /(2ir)dQ is the probability per unit of time 
of finding the photon in the infinitesimal spectral interval dO. We call 

Pxx (L,n) = ^\u(L,n)\ 2 (34) 

the spectral density of probability per unit of time of creating a photon-pair at the angular 
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frequencies uq ± |0|. Note that p xx is numerically equal to f x . 

When the condition (133]) holds, the pump scatters much less than one photon per mode. 
However, this condition is not strong enough to guaranty the validity of Eq. (1301) . For (I3"U|) 
being valid, the total probability Vt{L) = Tx J °° p xx (L, Q) dfl of scattering a photon during 
the time T must be much smaller than one. In order to find an analytical approximation of 
Vt{L), we neglect the term ^Pq x L in the argument of the sinc-function in Eq. (|32|) and find 
that Eq. fl30l) is an accurate approximation of the real quantum state of ligth if 



and the spectral width of the pump is Auj = 2tt/T, the square-root in Eq. ([3"5]) represents the 
number of independent modes in the fluorescence spectrum. Therefore, formula (I3l)j) shows 
that the expression (1301) is a good approximation of the quantum state of light only when the 
mean number of photons per mode is much smaller than the number of independent modes. 
However, photon-pair energy-entanglement over a wide spectral range is rarely desirable. 
If the Stokes and anti-Stokes photons are filtered, the condition for having two-particle 
energy-entanglement from FPS is that the mean number of photons per mode must be 
much smaller than the number of independent modes in the spectral interval allowed by the 
filtering process. 

D. Vector scattering 

When the pump wave is not polarized along an optical axis, not only scalar but also 
vectorial scattering can take place. The quantum state of light at the output of the fibre is 
then given by Eq. ffT6l) . The two-photon amplitudes £,ij(L, Q), with (i, j) 6 {x, y} 2 , are given 
in Eqs. ( ]29l) . Amplitudes £ xx and ^ yy correspond to a scalar scattering on the x and y axes, 
respectively. The remaining two amplitudes £ xy and £ yx correspond to a vectorial scattering 
in which two pump photons with orthogonal polarizations are annihilated. If we chose 
APi > (convention), the x-axis is the slow axis. Therefore (see Tab. [[]), £ xy corresponds 
to the scattering of the most energetic photon of the pair (the anti-Stokes photon) on the 
slow axis, while £ yx corresponds to the opposite process where the most energetic photon is 
polarized along the fast axis. 




(35) 



Since the spectral width of the Stokes (or anti-Stokes) spectrum is Af2 sca i/2 « a/27t/(|/3 2 |£) 
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As in the scalar case, the final state ^(L)) is an energy-entangled state. However, some 
polarization-entanglement is also present. The kind of polarization-entanglement that can 
be generated using HB fibers depends on the spectral shapes of the two-photon amplitudes, 
that in turn depend on the fiber parameters, the power of the pump, and its polarization 
state. We investigate the possible cases by first calculating the average spectral densities of 
photon flux, f x (L, fl) and f y (L, O), generated on the x and y axes, respectively. 



1. Spectral density of photon-flux 

Using Eq. (1311) (and a similar relation for the y-axis), one finds 

f x (L,n) = i_(|^(L,fi)| 2 + e(fi) |^(l,o)| 2 

+ |^(L,-0)| 2 ), (36a) 

f v (L,n) = ^% y (L,n)\ 2 + Q(n) \z yx (L,n)\ 2 

Z7T 

+ e(-n) |^(l,-o)| 2 ), (36b) 

where 9(0) = for Q < and 0(0) = 1 for O > 0. 

The functions f x (L,Q) and f y (L,Q) are plotted in Figs. [2](a) and^b), respectively, for 
typical parameters of HB silica fibers and an equipartition of the total pump power P 
between the optical axes (P 0x = P 0y = Po/2). Light gray and dark gray curves correspond 
to normal and anomalous dispersion; the other parameters are otherwise the same, including 
the absolute value of /3 2 - Three propagation lengths are considered. In all cases, L L n \ = 
min[L^j, L^], where L J nl = l/(7Poj) in the non-linearity length corresponding to the j-axis, so 
that the first-order perturbation formulas (|36l) are accurate approximations. An important 
feature of Fig. [2] is that the scalar and vector FPS take place in well separated frequency 
ranges. The broad central peak around Q = is due to scalar scattering, as described in 
Sec. IIVCI Two separate scalar FPS processes take place of the slow and fast axes. The 
narrow spectral peaks at higher frequency detuning 

lO^ctl w^(l-a) (37) 
IPs I 



are due to the vector scattering. In Eq. ([3 



& x 7-Po , Qfi s 

a = ~Kj3T- (38) 
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FIG. 2: Photon-flux spectral densities f x (L,£l) (a) and f y (L,Q) (b) as a function of $7 for prop- 
agation distances L = 100 m, L = 200 m, and L = 300 m. The other parameters used in this 
figure are 7 = 3 /(W km), P 0x = P 0y = 150 mW, Aft = 200 ps/km, j3 2 = ±15 ps 2 /km. Light 
gray and dark gray curves correspond to the normal (/?2 = +15 ps 2 /km) and anomalous dispersion 
(/?2 = — 15 ps 2 /km), respectively. 



It has been shown in 



351 ] that the dynamics of the fields in a HB fiber can be discussed in 



terms of this single dimensionless parameter. The situation displayed in Fig. [2] correspond 
to a < 1. As seen from the figure, changing the sign of the group-velocity dispersion 
has the same effect as interchanging Q and —Q. The narrow central peak at Q = is 
also due to vector scattering, but being positioned at the pump frequency, it is difficult to 
observe in practice. Only one of the amplitudes £ xy (L,Q) and £, yx (L, O) takes significant 
values for Q > [see the definitions fl29b ) and (I29H )]. In the normal dispersion regime, 



17 



N 
X 
H 

'w 

X 
13 



O 

h 

>> 

1 1 

tu 

I 

u 

CD 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 

L=45cm (L / L = 0.30) 



Slow axis 
■ Fast axis 




(L/L n| = 0.20) 



(L/L =0.10) 



W -30 -20 -10 10 20 30 
Angular frequency detuning Q [radTHz] 

FIG. 3: Photon-flux spectral densities f x (L,£l) (light gray) and f y (L,Q) (dark gray) as a function 
of 17 for propagation distances L = 15 cm, L = 30 cm, and L = 45 cm. The other parameters 
used in this figure are 7 = 36 /(W km), P 0x = P 0y = 20 W, Aft = 400 ps/km, (3 2 = -139 ps 2 /km 
(a = -1.25). 



vector FPS is due to non-zero values of £ yx in the neighborhood of the angular frequency 
(|37|) . In the anomalous dispersion regime, £ yx is zero for any Q > 0; £ xy takes significant 
values in the neighborhood of the angular frequency (1371) . instead. Tab. [T] explains the 
polarization properties of the scattered photons. The spectral width of the vector FPS peaks 
is Af2 vect w 4-71"/ (AfliL). As illustrated in Fig. [21 very monochromatic (Af2 vect w 0.2 radTHz) 
photon pairs can be generated using vectorial FPS. With picosecond pump pulses, Fourier- 
transform limited photon pairs could be generated without any external filtering scheme. 

As the value of the parameter \a\ increases, the frequency ranges of the scalar and vector 
scattering become closer and can even merge. In the Pq x = Pq v = Pq/2 case, this happens 



for |a| > 1. This behavior has been recently observed in a photonic-crystal fiber 35] in 
a modulation instability regime. Fig. [3] shows the spectral densities of photon-flux for a 



351 ]. but in a regime where L L n \ 



fiber with parameters similar to those in experiment 
min[L^j, L^J holds. Scalar and vector scattering occur in the same frequency range. One 
cannot distinguish photons from these different processes by their energy anymore. 
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2. Polarization- entanglement 



Energy-entanglement is a feature common to any, x or X^ 3 \ parametric photon-pair 
generation process. In contrast, intrinsic polarization-entanglement is only encountered in 
special circumstances. In x parametric down- conversion, for instance, it requires a type-II 



phase-matching [36] . In this paragraph, the kind of polarization-entanglement achievable 
using the FPS process in HB fibers is described. To separate polarization-entanglement 
from energy-entanglement, let us consider that the light coming out of the fiber is processed 
through a filtering apparatus that only selects, from the full spectrum, one pair of correlated 
modes with frequencies Uq ± Q and spectral width Au = 2tt/T (where T is, as before, the 
pump duration). We also consider, for simplicity, that P 0x = P 0y = Pq/2. 

Let's first consider the \a\ <C 1 case illustrated in Fig. [21 If Q corresponds to the center 
of the narrow vector scattering peaks, the quantum state of the generated photon pair is 

|lfcl-n> (39) 

for normal dispersion, and 

llSUn) (40) 
for anomalous dispersion. The generated photons always have opposite polarizations but 
they are not entangled. The polarization of the photons is correlated to their energy. 

Consider now that Q is somewhere in the frequency range corresponding to the scalar 
scattering. In that case, the generated photon pair will be in the quantum state 

^( • e* lol y „}). (41) 

where 

= 2(0oy-0o«) (42) 

The state (T4T1) is a maximally entangled two-particle state (Bell state). The relative phase 4> 
only depends on the polarization of the pump field. It can be controlled using a quarter-wave 
plate placed in the pump beam. Since the polarization of each photon of the pair can also be 
rotated individually after they exited the fiber, any Bell state can be produced. If the pump 
powers propagating on the optical axes are not equal (Pq x ^ Poy), the two components of 
the coherent superposition (j4Tj) do not have the same amplitude. In that case, the state is 
not maximally entangled. The fraction of the total pump power P that is polarized along 
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the slow axis can be controlled using a half-wave plate. Therefore, a full control of the pump 
polarization (quarter-wave and half-wave plates) permits to modify the relative phase and 
the strength of the entanglement. 

Consider now that the spectral range of scalar and vector scattering overlap. If the 
photons in the modes Uq ± Q can be produced by either the scalar or the vector process, the 
amount of entanglement will be reduced. In the extreme case depicted in Fig. [3] (\a\ > 1), 
the state 

(|l£lf. n >+e* |l^ n »+2x (|l£) + e<* |1&>) (|l^>+e«* |lf. n >) (43) 

is produced if Q is small enough (plateau in Fig. [3]). This state is partially entangled. 

Generating any possible polarization-entangled Bell states (T4T!) in a HB optical fiber 
is a very attractive prospect for quantum information processing. A first experimental 



demonstration has been realized very recently [17|, 118] (in a counter-propagating geometry 
to avoid any walk-off of the x- and y-polarized photons). The other processes described in 
this section have not been demonstrated yet. The process that producing the states fl39|) 
and (140]) is very interesting from an experimental point of view because it generates photon 
pairs that are very narrow-band (see Fig. [2]) and easy to separated (at the fiber end) using 
the polarization degree of freedom. 

V. LOW-BIREFRINGENCE FIBERS 

We now examine the peculiarities of FPS in fibers having a low birefringence (LB fibers). 
In that case, the group- velocity mismatch can be neglected (fii x = (3\ y = In contrast 
with the HB case, the last terms in the right-hand side of Eqs. ( 120]) cannot be dropped 
because the beat length 27r/A/3o may have the same order of magnitude as the other relevant 
length scale (the non-linearity length, for instance). The propagation equations are 

8A X dA x .p 2 d 2 A x . / 2 \ 

— + Pi-zr + l ^^r = % l A l A x + ~ A \K A x 



+ tlA 2 y Ale~ 2iA/3QZ , (44a) 



<):. " i)l ' ■> \ r 3 " 

,7 
3 

^y + B^ + i^i^y - iJtfA + 2 a^a\a 



+ i^-A 2 x A\e +2i ^ z . (44b) 
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As explained at the end of Sec. 1111} only FPS due to a monochromatic pump wave polarized 
along an optical axis will be studied here. We arbitrary choose this axis as the x-axis. The 
cartesian components of the envelope of the injected monochromatic pump are 

4,(0, t) = v 7 ^ e i6 °; A y (0, t) = 0. (45) 

Whether the x-axis is the slow or fast axis depends on the sign of A/3q. The monochromatic 
solution of Eqs. (j44H that is compatible with the initial condition (145 p is 

A x (z, t) = ^Po e ie ° e^ Poz ; A y (z, t) = 0. (46) 

As in the high-birefringence case, we want to linearize the propagation equations in the 
neighborhood of the monochromatic solution (|46j) . Therefore, we introduce the ansatz 

A x (z, t) = [y/F Q + u x {z, t)) e w ° e^ Poz , (47a) 
A y (z,t) = u y (z,t) e^ Poz (47b) 

and inject it in (jUJ. Linearizing the propagation equations with respect to u x and u y results 
in: 

du x du x .f3 2 d 2 u x . , n 

du y <9w y .I3 2 d 2 u y 

+l 1 Pq Uy e ^o ^(AA+hfl))*. (48b) 
3 

These equations show that the perturbations u x and it^ are uncoupled. Consequently, no 
correlation is expected between photons generated on the x- and y-axes. 



A. Coupled-mode equations 

Using Eqs. (j3J), (T4]), ffT9l and (jUj), the mode expansion of the fields u x and u y is obtained: 



x/ e— da ,49a) 

= y a .(^) e iB ^>- iQt dn, (49b) 
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with B(Q) = flifl + (/3 2 /2)0 2 . Injecting the expansions (|49|) in Eqs. fljgp , we obtaine the 
following coupled-mode equations: 

-^-a x (z,n) = i 7 P e i29 ° 
oz 

x e-ifr&+»r*>)* at(z,-n) (50a) 
^a y (z,Q) = z|P e^ 

x e -i(/3 2 ^-f 7 Po-2A/3o), fl t (Z) (50b) 

Note that the first equation is the same as equation (j27b ) for P 0y = 0. We know, from the 
study of Sec. IIV C\ that this corresponds to a scalar scattering on the x-axis, and will not 
repeat this analysis here. The interesting feature of FPS in LB fibers comes from Eq. (|50b ). 
which described a new kind of vectorial FPS, as can be seen from the structure of the 
infinitesimal generator: 



h 

2 



G lh (z) = '-^Po ( e 2ie ° I e -i(^ 2 +^Po)z a t( Z;fi ) a t( Zj _ fi ) dfi + H.c. 



4h Po ' e "° 



-i(f3 2 n 2 -hP -2Ap )z _t, 



aUz, n)aUz, -ft) dft + H.c. . (51) 



B. Four-photon scattering 

Using the infinitesimal generator (|5"T|) and moving to the interaction picture, the quantum 
state of light at the output of the fiber is found to be 



POO 

\^{L)) = |0)+ / £ ra (L,Q) |l£,l* n )dft 
Jo 

POO 

+ / £ yy (L,n) |i^i^ n )do, 

Jo 



(52) 



where £ xx is given by Eq. (I29a) (with P 0x = P Q ) and 



xsmc 



/? 2 ^ - - 7 P - 2A& } 



(53) 



The x-polarized pump not only scatters x-polarized photons through the scalar scattering 
process but also y-polarized photons through a vectorial scattering process. This vectorial 
process is different from those encountered in the study of FPS in HB fibers: here two 
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FIG. 4: Photon-flux spectral densities on the slow (gray lines) and fast (black lines) axes of the 
fiber as a function of Q. In the panel (a), the dispersion is normal and the pump in polarized along 
the slow axis of the fiber. In the panel (b), the situation is opposite: the dispersion in anomalous 
and the polarization of the pump wave is along the fast axis. The parameters used in this figure 
are 7 = 3 /(W km), Pn = 1 W, |A/3n| = 2 m , and | /3s | = 5 ps 2 /km. Three propagation distances 
(L = 50 m, L = 100 m, and L = 150 m) are considered. 



linearly co-polarized pump photons give birth to a signal-idler pair polarized orthogonally 
to them. 

The spectral density of photon flux associated with scattering in a low-birefringent fiber 
in shown in Fig. HI In Fig. H^, /?2 > and the pump is linearly polarized along the slow 
axis (A/3q > 0). The central peak is due to scalar scattering. The two narrow far-detuned 
peaks polarized along the fast axis are due to the vector FPS associated to the amplitude 
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(|53|) . Note that if the pump is polarized along the fast axis (A/3o < 0), an anomalous 
dispersion regime (/^ < 0) is needed to observe a similar scattering on the slow axis, as 
displayed in Fig. Hb. In both cases, the detuning and the spectral width (full width at the 
first-zeros of the sinc-function) of the generated photon pairs are |fi vec t,ib| ~ a/2A/3 //52 
and Af2 vectj ib ~ (2n / L) / ^2j3 2 Af3o, respectively. When A/3 x /3 2 < 0, FPS only occurs at 
frequency detunings close to £1 = 0. This case is not shown in Fig. HI 

Photon-pair generation using vector FPS in LB fibers has not been demonstrated yet. 
Vector FPS in LB fibers has however an important advantage over scalar FPS in terms of 
signal-to-noise ratio. It has been shown that, photon-pair correlations can be one order of 
magnitude higher than in the scalar FPS case for angular frequency detunings \Q\ between 30 
and 90 radTHz J^. Furthermore, photons generated by vector FPS are very narrow-band: 
as shown in Fig. HI Af2 vectj ib ~ 0.2 THz. By pumping with picosecond pulses, Fourier- 
transform limited photon pairs could be produced without filtering. Monochromaticity and 
low Raman noise are the two advantages offered by vectorial FPS in LB fibers. 



VI. CONNECTION WITH MODULATION INSTABILITIES 

There is a connection between FPS from a single pump beam and the phenomenon called 
"modulation instability" in nonlinear fiber optics [31]. In this work, a special care has been 
taken to formulate the perturbation theory of FPS in a way that makes this connection 
obvious. 

A modulation instability is a nonlinear phenomenon that allows an initially continuous 
monochromatic wave to become modulated in time due to the propagation in the fiber. 
This effect can be seen as the result of a degenerate four-wave mixing (FWM) process that 
is self phase-matched by a balanced between birefringence, dispersion and non-linearity. The 
pump develops exponentially growing symmetric sidebands about the frequencies ojq ± VL rn 
satisfying phase- matching conditions. In the time domain, the envelope of the pump is 
therefore modulated at angular frequency fl m . 

The dynamics of modulation instabilities can be understood by solving the equation 
systems ( )27l) or (J50l) exactly. This method is used in classical nonlinear optics where ctj (z, Q) 
and a^(z,Q) are treated as classical spectral amplitudes. The analysis provides the phase- 
matching conditions, the frequency range of the instable (i.e. exponentially growing) modes 
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experiencing the parametric gain, the value of that gain, and the polarization of the growing 
sidebands. The development of the instability can be induced by a additional coherent probe 
signal at u ± Q m or even some incoherent optical noise travelling with the pump. If care is 
taken to eliminate these instability sources, modulation instability can still develop because 
the photon pairs produced by the FPS phenomenon populate the unstable modes. In that 
case, the modulation instability is called spontaneous, and has a purely quantum origin. 

It is important to note that although modulation instabilities can be trigger by the FPS 
process, their spectral properties are different. Let us illustrate that point in the case of a 
scalar modulation instability. The analysis of Eqs. ( 1271) with P Qx = P and Po y = shows 
that instability is only observed in anomalous dispersion regime (fli < 0). (In reality, an 
instability can also be observed in the normal dispersion regime if the pump wavelength 
is close to the zero- dispersion wavelength; this results from higher order dispersion effects 
that have been omitted here by limiting the Taylor expansion (JTTj) to second order.) The 



parametric gain g(£l) = a/(7-^o) 2 — (27P0 — |/^2|^ 2 ) 2 /4 only exists in the angular frequency 
range ]— 2^ r yP /\^\, 2y / 7P /|/^2|[- Its maximum value g max = 7P0 is reached for the angular 
frequencies O max = ±^2'yPo/\f32\ corresponding to a perfect phase- matching with the pump. 
For g{Q)L ^> 1, the photon flux that results from a spontaneous scalar modulation instability 
is well approximated by 

f x (L,0) = — -^-e 2 ^ L . (54) 
J v ' ; 2tt4 5 2 (Q) v ; 

Because of the exponential growth, the spectrum associated with modulation instability ex- 
hibits sharp peaks at detunings ±Jl max . This strongly contrasts with the broad fluorescence 
spectrum due to scalar FPS. The build-up of the modulation instability can be understood 
by comparing the spectral intervals over which parametric amplification and FPS take place. 
For a given pump power Pq, the width of the spectral interval in which scalar modulation 



instability develops is Af2 SM i = Aiy^P /\p 2 \- It is independent of the propagation length. 
In contrast, FPS takes place in a spectral interval Afi scal 2^/27r/(|/3 2 |-^), that depends on 
the propagation length L but is independent of the pump power to a good approximation 
(see Sec. IIV Cp . The ratio of these bandwidth is r = Af2 SMI /Af2 sca i « 0.8 x ^^P L. This 
shows that, in the beginning of the propagation {^PqL 1), the FPS spectrum is much 
broader that the amplification band. As light propagates further, the number of photons per 
mode increases and the FPS spectrum becomes narrower. When ^PqL rs 1 (the propagation 
length is equal to the non-linearity length L nl = I/7P0)) the width of the FPS spectrum is 
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FIG. 5: Graphical representation of the mechanism leading to the growth of the number of photon 
pairs in the field in the case of a scalar scattering. Each arrow represents a scattering event, i.e. 
the creation (right-pointing arrow) or the annihilation (left-pointing arrow) of a signal/idler pair. 
Column's labels represent the number of photon pairs in the field. Rows' labels indicate the order 
in the Dyson's perturbation series accounting for the scattering events represented by arrows in 
that row. 

almost equal to the width of the parametric amplification band. At that moment, each mode 
in the amplification band has been populated with about one photon per mode. Stimulated 
FWM can therefore take place and modulation instability develops. 

The discussion above is qualitative because, when the mean number of photons per mode 
approaches one, the scattering theory based on the first-order perturbation approach (FPS) 
is not accurate anymore. To describe what happens when the number of photons scattered 
from the pump is higher than one, higher-order perturbation theories based on Eq. (floT) with 
n > 1 can be used. Fig. [5] shows how the different perturbation orders (fl~5l) in Dyson's series 
contribute to the quantum state of light in the case of scalar scattering. 

The rows represent successive perturbation orders, while the columns represent the num- 
ber of photon pairs in the field. An arrow connecting a node with n photon pairs at per- 
turbation order k to a node with I photon pairs at perturbation order k + 1 means that 
the amplitude of the Z-pair state (at order k + 1) depends on the amplitude of the n-pair 
state (at order k) . The figure shows that limiting the Dyson's series to the k th order consists 
in approximating the quantum state of light by a state containing no more that k photon 
pairs. This approximation is therefore precise if the probability of having k photon pairs in 
the field is small. It also shows that the amplitude of a n-pair state at a given perturbation 
order k only depends on the amplitudes of the states with n — 1 and n + 1 photon pairs at 
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order k — 1. For instance, in the second order approximation, the quantum state is 

POD 

M = [l-o(L)]|0)+ / \l x n A-n) dft 

Jo 



+Uj ° O U(L,n) |l£,l* n )dn) . (55) 

The term in the second line of (|55p represents the possibility of emitting two photon pairs, 
and a(L) is the second-order correction to the amplitude of the initial "no photon-pair" 
state |0). The probability of creating at least one photon pair is equal to 2Re[a(L)]. This 
number only takes a finite value if the pump wave has a finite duration T. For ^Pq x L <C 1, 
2Re[a(L)] ~ Vt(L) [see Eq. (I35j) ] . One can easily compute that the mean number of photons 
generated during the time T in a spectral mode around the frequency detuning ft is 

n x (L, ft) = \U(L,n)\ 2 + \U{L,tt)\ 2 V T {L) + \£ xx (L,n)\\ (56) 

The first term corresponds to the first-order approximation while the second and third terms 
are second-order corrections. The second term describes a second independent FPS event. 
The probability of scattering a second photon in the mode ft is equal to the probability of 
scattering a first photon (regardless the mode) times the probability of scattering a photon 
in the mode ft. The third term is interpreted as a stimulated emission contribution. An 
additional photon in the ft mode is more likely if the first one as already been emitted in 
that mode. 

Higher-order perturbation theory provides important insights to the physics of photon- 
pair generation in fibers and the built-up of spontaneous modulation instabilities. However, 
it is clearly not an efficient way to compute the mean number of photons in a given spectral 
mode. In the case of scalar scattering, the quantum equations of motion (j2"Tj) . with P 0x = 
Pq and Pay = 0, can be solved analytically without any approximation. This solution is 
equivalent to summing the Dyson's series. The average photon-flux spectral density derived 
by that method is 

f x (L,Q) = |sinh(A(ft)L)| 2 , (57) 

where A (ft) = V(7 p o) 2 - (27P0 + /3 2 ft 2 ) 2 /4. In the limit -fP L < 1, the formula (T57|) 
reduces to ( |32l) . In the anomalous dispersion regime (fa < 0) and the ■yPoL ^> 1 limit, it 
gives Eq. (15^1) ; note that A (ft) can only be identified with the gain g (ft) when A (ft) G K. 
Eq. (1571) is a well known result for a parametric amplifier 28] . In the case of vectorial 



scattering, the quantum equations of motion, Eqs. ff27l) in the HB case or Eqs. f )50l) in the 
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LB case, can be solved without approximations only when only one scattering process con- 
tributes to the population of the signal (— fl) and idler (il) modes under consideration. For 
these cases, formulas similar to (1571) can then be derived. In HB fibers, different processes 
can scatter pump photons to the same spectral interval. As seen in Sec. IIVD[ this happens 
when the parameter \a\ [Eq. (1381) ] is higher than one. In that case, Eqs. (12 7p cannot be solved 
analytically. Perturbation theory, however, permits to study these cases. From a general 
point of view, solving the equations of motion exactly is not a good theoretical approach to 
the photon-pair generation problem because the physical process of FPS cannot be straight- 
forwardly derived from the mathematical solution. First-order perturbation theory, on the 
other hand, clearly exhibits the physical process of FPS, applies to any physical situation, 
gives accurate results for L <C L n \ and can be generalized (higher-order perturbation) the 
take into account multiple scattering events. 



VII. CONCLUSION 



FPS in nonlinear waveguides is an important physical process that enables photon-pair 
generation in well defined guided modes, with high rate and low noise. Most of the exper- 
iments to date used the scalar FPS process in which the pump photons and the scattered 
photons have the same polarization. In birefringent waveguides, vectorial FPS processes 
are also allowed: these vectorial scattering processes involve photons with different polar- 
izations. In this article, the theory of FPS in nonlinear, birefringent, and dispersive fibers 
is developed in the framework of the quantum theory of light. It gives the general method 
for studying different vectorial FPS processes, based on the first-order perturbation theory 
of scattering. 

Different FPS occur in high- and low-birefringence fibers. These cases have been stud- 
ied separately. In the high-birefringence case, photons with orthogonal polarizations can 
be generated. According to the process used for their generation, they can be either en- 
tangled in polarization or have their polarization correlated to their wavelength. In the 
low-birefringence case, photon pairs with polarization orthogonal to the pump field can be 
produced. This generation methods has the advantage that polarization can be use to sepa- 
rate the photons from the pump and that Raman noise is lower, as has been proved recently 
29l . l30fl . Our main concern has been the spectral properties of the generated pairs. We 
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showed that very monochromatic photon pairs can be generated using the vectorial FPS 
processes. Consequently, when picosecond pump pulses are used, it is possible to generated 
photon pairs with Fourier-transform limited spectra. We also noted that scalar and vector 
processes may sometimes scatter photons into the same spectral bands. 

In order to apply the theory to practical design of photon-pair sources, it may be desirable 
to extend it in different directions. For instance, the theory can be extended to encompass 
more complex dispersion relations, as required when working close to the zero-dispersion 
wavelength or/and with photonic crystal fibers. This can be done by modifying the dis- 
persion operator in the quantum nonlinear Schrodinger equations ( 120]) . Exotic dispersion 
can dramatically modify the spectral properties of the generated photons. Including the 
spontaneous Raman effect is important for estimating the signal-to-noise ratio of a source. 
This can be done by generalizing the nonlinear terms of the quantum nonlinear Schrodinger 
equations (1201) as explained in [30j. For deriving the quantum state of light at the output of 
the fiber and the mean photons fluxes, the method explained in this paper applies and can 
by followed step by step. For dealing with finite pump duration (pump pulses), the easiest 
method is to use discrete modes, as we did it this work. If the precise pump shape must be 
taking into account, the theory has to be modified more seriously. This can be done in the 
same way as in the scalar scattering case [19J], but that method is not rigorous. A better 
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method for dealing with pump pulses is to integrate ([20]) numerically, as explained in 

Understanding the physics of FPS in birefringent media is important for the present 
development of the field of quantum photonics. This work contributes to this rapidly growing 
field by presenting the principles underlying vectorial FPS and, as explained above, can be 
extended in many different ways to match the need of particular applications. 
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